Chaos in Symmetric Phase Oscillator Networks 
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Phase-coupled oscillators serve as paradigmatic models of networks of weakly interacting oscillatory units in 
physics and biology. The order parameter which quantifies synchronization was so far found to be chaotic only 
in systems with inhomogeneities. Here we show that even symmetric systems of identical oscillators may not 
only exhibit chaotic dynamics, but also chaotically fluctuating order parameters. Our findings imply that neither 
inhomogeneities nor amplitude variations are necessary to obtain chaos, i.e., nonlinear interactions of phases 
give rise to the necessary instabilities. 

PACS numbers: 05.45.Xt, 02.20.-a 



Introduction — Models of coupled oscillators describe var- 
ious collective phenomena in natural and artificial systems, in- 
cluding the synchronization of flashing fireflies, or supercon- 
ducting Josephson junctions, oscillatory neural activity and 
oscillations in chemical reaction kinetics |1|. In particular, 
phase-coupled models serve as paradigmatic approximations 
for many weakly coupled limit cycle oscillators \2. 3|. The 
Kuramoto model (and its extensions) provides the gold stan- 
dard in this field because it suitably describes the dynamics of 
a variety of real systems, is extensively studied numerically 
and reasonably understood analytically ||4jt5|. Each oscillator 
k with phase V'fc(i) € M/27rZ =: T on the 1 -torus T changes 
with time t according to 

for all fc € {1, . . . , A^}. For the original Kuramoto model 
the coupling function g has a single Fourier mode, g = sin. 
The dimension of such systems can be reduced to low di- 
mensions |6-8|, implying dynamics that is either periodic 
or quasi-periodic. For coupling functions with two or more 
Fourier components the collective dynamics may be much 
more complicated. For example, stable heteroclinic switching 
may emerge lUKTO). More irregular, chaotic dynamics of sys- 
tem ([TJ is observed for non-identical oscillators only lfmfT2l . 
raising the question whether inhomogeneities are necessary 
for the occurrence of such dynamics. To the best of our knowl- 
edge, the only hint that chaotic dynamics might exist for sym- 
metric phase-oscillator networks are attractors with irregular 
structure in phase space found recently in a system of = 5 
oscillators f9\. 

The complex order parameter 

1 ^ 

where i = \/— T constitutes an important characteristic for 
coupled oscillator systems. In particular, its absolute value 



\R{t)\ quantifies their degree of synchrony with \B.{t)\ = 1 
if all oscillators are in phase. For the original Kuramoto sys- 
tem the full complex order parameter (|2]) acts as a mean field 
variable enabling closed-form analysis lfT3l . 

In homogeneous, globally coupled systems ([T} it remains 
unknown whether or not there exists any coupling function g 
that gives rise to chaotic order parameter fluctuations. Syn- 
chronous solutions, anti-synchronous splay states as well as 
the dynamics of cluster states (where ipk — (fij for at least 
one k j) all yield a periodic complex order parameter R{t). 
For any invariant solutions on tori, R{t) is either periodic or 
quasiperiodic. The most irregular dynamics of R{t) observed 
so far is due to heteroclinic cycles where R{t) is non-periodic 
as it 'slows down' each cycle. Even if chaotic dynamics does 
emerge within the system, it may average out due to sym- 
metry, possibly resulting in a regular dynamics of the order 
parameter 

In this Letter, we answer the question whether chaotic dy- 
namics, and moreover, chaotic order parameter fluctuations 
may arise for some g even in the absence of inhomogeneities 
in ([T]). We show that indeed chaos is not possible for A^ < 4. 
By contrast, for A^ = 4, chaotic attractors can appear in a spe- 
cific family of coupling functions g. Interestingly, attractors 
of all theoretically possible symmetries exist and we provide 
further examples of attracting chaos for N ~ 5 and N ^ 7. 
The existence of chaos for infinite families of A^ > 4 implies 
that chaos occurs in systems with certain N > Nq for any 
Aq e N and suggests that chaos is likely to occur in many 
high-dimensional systems with a suitable choice of g. 

No chaos for N ^ 2 and N ^ 3 — Let be the A^- 
dimensional torus and let Sn denote the group of permuta- 
tions of A^ symbols. Suppose M is a differentiable manifold 
and let F be a group that acts on M. Recall that a vector field 
X on M is called F-equivariant if X 'commutes' with the ac- 
tion of r, i.e., X o 7 = 7 o X for all 7 e F where 7 denotes 
the induced action on the tangent space. 

Equivariance implies restrictions of the dynamics specified 
by the vector field. We study the dynamical system on 
given by the ordinary differential equations ([T]). Let us hence- 
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Figure 1 : Trajectories of the absolute value of the order parameter 
\R{t)\ fluctuate chaotically (A'^ — 4). Two trajectories with small dif- 
ference in initial condition diverge from each other (coupling func- 
tion ^ for parameter values rji — 0.1104 and rj2 — 0.5586). 



forth assume that the system is homogeneous, i.e., lo^ = lj for 
all fc S {1, . . . , N}. This systemis Sm'xT^ equivariant where 
Sm acts by permuting indices and through a phase shift. 
Recall some basic properties of this system Q. Introducing 
phase differences ipj :— (pj—(pi for all j S {1, . . . , N} elimi- 
nates the phase-shift symmetry. Write (p = {(pi, . . . , (pjv) and 
4> = (i/'ij • ■ ■ J V'at)- The reduced system on T^^^ is given by 




■0fc) 



N 

E 

k=l 



gi-iJ-'k) 



(3) 



forallje{2,...,iV}. 

For any partition P = {Pi .... , f ,„} of {1, .... iV} (that is 

Pr C {1, . . .,N}, \J]1, P, ={!,.. .,N), and n = 
for r ^ s) the subspaces 



Fp := {ip \ j,k G Pr for any ; 



= (4) 



are flow-invariant. The subspaces divide ^ in {N — 1)! 
invariant {N — 1) -dimensional simplices |2|; one of which 



C ;= {i/. I = V-i < ^-2 < ■ ■ ■ < < 27r} C T 



(5) 



we refer to as the canonical invariant region. There is a 
:= Z/NZ symmetry on the canonical invariant region and 
the 'splay state' (the phase-locked state with t/jj = 2TTj/N) is 
the only fixed point of this action at the centroid of this region. 

The reduction of symmetry has implications for the exis- 
tence of chaos in low dimensions. For N = 2 and = 3 
the phase space of the reduced system is a one, resp. two- 
dimensional torus. This means that by the Poincare-Bendixon 
theorem lfT4l chaos is not possible in these systems for < 4. 

Chaos and symmetry for N — A — We choose a 
parametrization of the coupling function g in ([TJ by consid- 
ering a truncated Fourier series 



div) = ^o,k cos{k(f + £,k) 



(6) 



k=l 



In particular, we restrict ourselves to the two parameter family 
given by the parametrization 

(6,6,6,^4) = {vi,-vi,m +V2,vi + V2) (7) 
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Figure 2: Chaos for A'^ = 4. Maximal Lypunov exponent is positive 
in a region of parameter space. The initial condition was fixed and 
the coupling function parametrized by 



where r/i and 772 are real valued parameters and oi = 
—2, 02 = —2, 03 = —1, and 04 = —0.88 are constants. 

For = 4, chaotic attractors do indeed exist. The dy- 
namics of the absolute value of the order parameter exhibits 
chaotic fluctuations and exponential divergence of trajecto- 
ries, cf. Figure [T] To explore parameter space, we calculated 
the maximal Lyapunov exponent Amax from the variational 
equations ll20l . There are regions in (771, 772) -parameter space 
in which Amax is greater than zero, cf. Figure |2] As might 
be expected, there is fine structure in this region, for example 
islands where the trajectory converges to a stable limit cycle. 
Lines of period doubling cascades |21 1 bound the chaotic re- 
gion and end in a homoclinic flip bifurcation with an inclina- 
tion flip 1 15 1 (details not shown). Exploring initial conditions 
revealed the coexistence of chaotic attractors and stable limit 
cycles in part of the chaotic region. 

Chaotic attractors in equivariant dynamical system can ex- 
hibit symmetries themselves. Let ^ be a chaotic attractor as 
defined in [161 , i.e., a Lyapunov-stable, closed, and connected 
set that is the w-limit set of a trajectory, for a dynamical sys- 
tem on a manifold M given by a F-equivariant vector field. 
The subgroup Stab(A) := {7 e F | 7(a) = a for all a € A} 
is the group of instantaneous symmetries of the attractor, 
i.e., at any point in time the action of Stab(^) keeps ev- 
ery point in A fixed. Furthermore, we define := 
{7 e r I 7(A) = A} to be the set of symmetries on average, 
and we have Stab(A) C as a subgroup. 

The subdivision of the phase space by flow-invariant re- 
gions restricts the possible symmetries of chaotic attractors. 
The possible symmetries on average of any chaotic attractor 
A C C with trivial instantaneous symmetries Stab(A) = {1} 
are limited to subgroups of Zat since they are contained one of 
the invariant simplices (C or one of its images under the group 
action) with that symmetry. For N — 4, any chaotic attractor 
of this type must have trivial instantaneous symmetry. Thus, 
the possible symmetries on average are limited to subgroups 
of Z4, i.e., C Z4. In fact, we have found examples 
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(a) S(A) = 1 (b) S(A) = (c) S(yl) = Z4 



Figure 3: All possible symmetries of the chaotic attractors for different parameters for A'^ = 4. We have rji = 0.138 in panel (a), rji — 0.0598 
in panel (b), rji — 0.1104 in panel (c), and 772 = 0.5586 in all panels. The projection is a F = S4 equivariant map, xi — sm((p2 — (^4), 2:2 = 
sin((^i — (/P3), and X3 = \R\. 



of chaotic attractors for each possible symmetry in systems 
of iV = 4 and coupUng functions given by (|7]) (Figure [3]l. 
Note that this definition of attractor is somewhat restrictive — 
Milnor attractors may display a wider range of symmetries in- 
cluding different instantaneous symmetries at the same time. 

Chaos for N > 4 — Analyzing the same region of pa- 
rameter space for > 4 yields attracting chaos in systems 
of A^ = 5 and N — 7 oscillators in large regions. Figure]?] 
shows an overlay of regions for three different A^; regions are 
shaded where the Lyapunov exponent exceeds 0.01 and darker 
areas indicate that several A^ satisfy this condition. Clearly, 
there is a single coupling function for which attracting chaos 
is present for all A^ = 4, A^ = 5 and, N = 7. Intriguingly, we 
did not find chaotic attractors for any A^ g {6, 8, 9, . . . , 13} in 
the entire region of parameter space considered in Fig. |4] 

The parametrization of the coupling function by a truncated 
Fourier series raises the question how many Fourier compo- 
nents the coupling function needs to contain for chaos to oc- 
cur For A^ = 5 we also measured positive Lyapunov expo- 
nents when the coupling was chosen to be through the simpler 
coupling function 5 ((ys) = — 0.2 cos((/3 + 771) — 0.04cos(2<y9 — 




Figure 4: Overlapping regions where the maximal Lyapunov expo- 
nent is greater than 0.01 for A'^ G {4, 5, 7}. The darker the color, the 
more for which the condition holds. For parameter values around 
(0.115, 0.06) there is a region where there is chaos for all these three 
N. 



772) as in ||9|. Hence, in dimension five, coupling functions 
with only two Fourier components suffice to generate chaotic 
dynamics whereas for A^ = 4, we did not find an example 
with less than four components. 

From the above, it is clear that for systems of size A^ = 
KM with K £ {4, 5, 7} there are chaotic invariant sets ly- 
ing in flow-invariant subspaces for coupling functions yield- 
ing positive Amax, cf. Figure ]4] For instance, for K = A, 
these spaces are given by partitions P — {Pi, . . . , P4} with 
\Pj\ = M for j e {1, . . . , 4}. For A^ large we similarly 
calculated a positive maximal Lyapunov exponents for the 
system reduced to asymmetric 4-cluster states given by par- 
titions P = {Pi, . . . , F4} with \Pi\/N = 1/A + q and 
|P,|/A^ = l/4-(7/3forj e {2,3,4} as depicted in Figure]!] 
However, these chaotic invariant sets in invariant subspaces 
close to the symmetric cluster state may be transversally re- 
pelling, possibly yielding non-chaotic long-term dynamics. 

Discussion — Inhomogeneities or asymmetries are thus not 
necessary for collective chaotic dynamics to appear in system 
([TJ, and even chaotic order parameter fluctuations emerge in 
the presence of full Sn x T-symmetry. We highlighted that 
for certain coupling functions chaotic attractors exist for sev- 
eral A^. However, the regions in parameter space for which 
chaotic attractors exist vary drastically, cf. Figure ]4] For cer- 
tain coupling functions there are chaotic invariant sets lying 
in flow-invariant subspaces that correspond to the symmet- 



0.15 




-0.005 -0.0025 0.0025 0.005 

Q 



Figure 5: Positive maximal Lyapunov exponents for asymmetric 4- 
cluster states for large systems, = KM + L{q, M) ^ 4. Here 
q parametrizes the deviation from the symmetric cluster state and L 
is the corresponding integer dimension (coupling function (J?} with 
771 = 0.1104 and r?2 = 0.5586). 
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lie and near-symmetric cluster states but these may not be 
transversally attracting. The question remains whether there 
are couphng functions giving rise to chaotic sets that are ac- 
tually attracting for = 6 and N > 8. Moreover, is there a 
'universal chaos function' in the sense that there is a coupling 
function for which there is some Nq E N such that there ex- 
ists a chaotic attractors for all (or at least an infinite number 
of)N > NqI 

For coupling functions with only one Fourier component, 
finite dimensional systems and the continuum limit are re- 
lated; the dynamics for both finite N and in the continuum 
limit reduces to effectively two-dimensional dynamics 171 El 
preventing the occurrence of chaotic trajectories. Is chaos pos- 
sible in the continuum limit for more complicated coupling 
functions? If so, how would such a result relate to chaos in 
the finite-dimensional systems we have studied here? 

For finite systems, attracting chaos in the system does not 
necessarily imply chaotic dynamics of the order parameter 
since there could be chaotic fluctuations for example in an 
'antiphase state.' The converse, however, holds. Addition- 
ally, observed chaotic fluctuations of the order parameter can- 
not necessarily be traced back to inhomogeneity in the system 
(possibly through an experimental setup, cf. [17]) because, as 
shown above, even fully symmetric systems can support such 
dynamics. When considering the continuum limit, the prob- 
lem of these implications becomes more subtle and will re- 
quire further investigation. 

The coupling function we considered above are written in 
terms of a truncated Fourier series. As discussed above, the 
number of Fourier components is relevant for the dynamics. 
An alternative approach would be to consider suitable piece- 
wise linear functions. For iV = 4, we find that systems with 
piecewise linear g also exhibit positive maximal Lyapunov ex- 
ponents (not shown). Finding a suitable basis for the space of 
coupling functions might be a way to explain some of the dy- 
namical features that were observed. 

Coupled phase oscillators are a limit of weakly coupled 
limit cycle oscillators [21. In globally coupled identical 
Ginzburg-Landau oscillator ensembles, chaotic dynamics can 
be observed ifTSlfTOl . However, it was thought that the ampli- 
tude degree of freedom is crucial for the emergence of such 
dynamics. Our results show that this is not the case and, 
moreover, suggest that chaotic mean field oscillations are also 
present in a large class of higher-dimensional symmetrically 
coupled limit cycle oscillators with a rich possible range of 



chaotic dynamics. 
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